Astronomy & Astrophysics manuscript no. 14767 


©ESO 2010 


June 3, 2010 





Multi-frequency study of Local Group Supernova Remnants 

The curious case of the Large Magellanic Cloud SNR J0528-6714 

EJ. Crawford^ M.D. Filipovic^'^, F. Haberl^, W. Pietsch^, J.L. Payne\ and A.Y. De Horta' 

' University of Western Sydney, Locked Bag 1797, Penrith South DC, NSW 1797, Australia 

e-mail: e.crawford@uws.edu.au, m.filipovic@uws.edu.au, a.dehorta@uws.edu.au 
~ Max-Planck-Institut fiir extraterrestrische Physik, Giessenbachstrasse, D-85741 Garching, Germany 
e-mail: fwh@mpe.mpg.de, wnp@mpe.mpg.de 

Preprint online version: June 3, 2010 

^ , ABSTRACT 

Ch ' Context. The SNRs known in the Local Group show a variety of morphological structures that are relatively uncorrelated in the 

^ I different wavelength bands. This variety is probably caused by the different conditions in the surrounding medium with which the 

' ' , remnant interacts. 

' Aims. Recent ATCA, XMM-Newton and MCELS observations of the Magellanic Clouds (MCs) cover a number of new and known 

SNRs which are poorly studied, such as SNR J0528-6714. This particular SNR exhibits luminous radio-continuum emission, but is 
one of the unusual and rare cases without detectable optical and very faint X-ray emission (initially detected by ROSAT and listed 
as object [HP99]498). We used new multi-frequency radio-continuum surveys and new optical observations at Ha, [Sii] and [Oiii] 

O wavelengths, in combination with XMM-Newton X-ray data, to investigate the SNR properties and to search for a physical explanation 
' for the unusual appearance of this SNR. 

(~| I Methods. We analysed the X-ray and Radio-Continuum spectra and present multi-wavelength morphological studies of this SNR. 

Results. We present the results of new moderate resolution ATCA observations of SNR 10528-6714. We found that this object is a 
I ■ typical older SNR with a radio spectral index of (!r=-0.36±0.09 and a diameter of D=52.4± 1 .0 pc. Regions of moderate and somewhat 

O ' irregular polarisation were detected which are also indicative of an older SNR. Using a non-equilibrium ionisation collisional plasma 

^ I model to describe the X-ray spectrum, we find temperatures kT of 0.26 keV for the remnant. The low temperature, low surface 

C/2 , brightness, and large extent of the remnant all indicate a relatively advanced age. The near circular morphology indicates a Type la 

■ event. 

Conclusions. Our study revealed one of the most unusual cases of SNRs in the Local Group of galaxies - a luminous radio SNR 
without optical counterpart and, at the same time, very faint X-ray emission. While it is not unusual to not detect an SNR in the 
^ , optical, the combination of faint X-ray and no optical detection makes this SNR very unique. 

, Key words. Galaxies: Magellanic Clouds ~ ISM: supernova remnants ~ individual: SNR J0528-6714 
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• 1. Introduction Here, we report on multi-frequency observations of a 

previously known and intriguing LMC supernova remnant. 

O The Large Magellanic Cloud (LMC), with it s low for eground ab- SNR J0528-6714 was initially suggested as a candidate by 

rZ^ sorption and relative proximity of -50 kpc (Idi Bene detto 200^, & Mdl^ im3- | Mathewson et alj 11985) classified it as 

• • offers the ideal laboratory for the detailed study of a complete ™ ^NR based on observations made with the Molonglo Radio 

> sample of objects such as supernova remnants (SNRs). The Observatory, at 36cni (v=843MHz), and they noted no optical 

proximity enables detailed spatial studies of the remnants, and counterpart. Filip ovic et al. | ^1998) added further confirmation, 

the accurately known distance allows for analysis of the energet- ^ith a set of radio-continuum observations (with the Parkes tele- 



X 



■ ics of each remnant. In addition, the wealth of wide-field multi- scope) over a wide frequency range and estimated a steep spec- 

wavelength data available, from radio maps to optical emission- t™m with powerlaw index a = -0.79 ± 0.15.|BlairetalJ (12Q06D 

line images and broad-band photometry to global X-ray mo- observed SNR J0528-6714 but reported no detection at far- 

saics, provides information about the contexts and environments ultraviolet (FUV) wavelengths. | Haberl & Pietsc h a999) associ- 

in which these remnants ai-e born and evolve. ^ted an X-ray countei-part to this radio SNR from ROSAT PSPC 

observations. According to the entry number in their catalogue. 

It is possible to obtain a relatively complete sample of SNRs the ROSAT source is named [HP99] 498. 
in the LMC and not only study the global properties of the sam- 
ple but also study the subclasses in detail (e.g., sorted by X-ray 
and radio morphology, diameter or by type of the supernova 

progenitor). Toward this goal, we have been studying SNRs in 2.1. The ATCA radio-continuum observations 
the Magellanic Clouds (MCs) in greater detail using combined 

optical, radio, and X-ray observations. Today we know over We observed SNR J0528-6714 with the Australia Telescope 

40 confirmed SNR s in the LMC and another 35-40 candidates Compact Array (ATCA) on 6* April 1997, with an array con- 

(iPavne et alJl2008h . figuration EW375, at wavelengths of 6 and 3 cm (v=4790 and 



2. Observational Data 
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8640 MHz). The observations were done in so called "snap- 
shot" mode, totaling ~1 hour of integration over a 12 hour pe- 
riod. Source 1934-638 was used for primary calibration and 
source 0530-727 was u sed for seco ndary ca libration. The miriad 
(ISault & Killeenll2006h and KARMA (lGoochl2006.) software pack- 
ages were used for data reduction and analysis. More infor- 
mation about the observing procedure and other sources ob- 
served during this session can be found in Bojicic et al. (2007|), 
ICrawfordetall (l2008allbh and lCaiko et all (l2009h . 

Baselines formed with the 6* ATCA antenna were excluded, 
as the other five antennas were arranged in a compact config- 
uration. The 6 cm image (Fig. [TJ has a resolution of 39"x31" 
at position angle 50° and the r.m.s noise is estimated to be 
0.4 mJy/beam. Due to the signal to noise restrictions and the 
size of the remnant, no reliable image could be prepared at 3 cm. 

We also used all available radio-continuum images of the 
LMC. These are composed of observations at several radio fre- 
quencies having moderate r esolution at 36 cm (v=843 MHz, 
MOST; [T urtle & Amvlll99lh and 20 cm (v=1400 MHz, ATCA; 
iHughes et alJi2007ir "^ 

2.2. XMM-Newton observations and data reduction 

XMM-Newton serendipitously observed SNR J0528-6714 on 
28* January 2007, for a total of about 20 ks (observation ID 
0402000601) at an off-axis angle of ~6'. The observation was 
perfor med with the EPIC instruments (PN ari d two MOS cam- 
eras, IStriider et alj 120011; [Turner et"al] 1200 Ih in imaging read 
out mode. Thin optical blocking filters were used to optimise 
observations of the tar get, the superso ft X-ray source candi- 
date RX J0529.4-67 1 3 iKahabka et al.1 l2008l) which is associ- 
ated with the planetary nebula SMPL69. The data were anal- 
ysed using the analysis package XMMSAS version 8.0.0. After 
removal of intervals with high background activity, we obtained 
a net exposure time of 17.3 ks (EPIC-PN). 

For morphology studies in comparison with radio-continuum 
and optical images XMM-Newton EPIC images were pro- 
duced in the standard energy bands 0.2-0.5 keV, 0.5-1.0 keV, 
1.0-2.0 keV and 2.0-4.5 keV. The images show faint diffuse 
emission at the position of the SNR and a point-like source near 
its centre at RA(J2000)=5''28'"17.9', DEC(J2000)=-67°14'01" 
(statistical error 1", systematic uncertainty 2-3")- 

To investigate if the point-like source is related to the SNR 
we extracted spectra for both sources separately. For the point- 
like source we used a circular extraction region with 20" ra- 
dius. For the SNR the outer radius of the extraction region (cen- 
tre at RA(J2000)=5''28'"19.8', DEC(J2000)=-67°14'21") was 
100", excluding a circle around the point-like source with ra- 
dius 25". EPIC X-ray spectra were extracted for PN (single + 
double pixel events, corresponding to a PATTERN 0-4 selec- 
tion) and MOS (PATTERN 0-12), excluding bad CCD pixels 
and columns (FLAG 0). We used XSPECQ version 12.5.1 for 
spectral modelling. The three EPIC spectra from the SNR were 
fitted simultaneously, allowing only a renormalisation factor to 
account for cross-calibration uncertainties between the instru- 
ments and area losses for the MOS extraction as the region cov- 
ered by the SNR is near the edge of the inner MOS CCD in 
both cameras. To account for photo-electric absorption by inter- 
stellar gas, two hydrogen column densities were used. The first 
represents the foreground absorption in the Milk y Way, fixed 
at 6x 10^" cm"^ assuming elemental abundances of IWilms et al.l 
(12000.) . The second considers the absorption in the LMC (with 

' http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 



Table 1. Integrated Flux Density of SNR J0528-6714. 





Frequency (MHz) 
Wavelength (cm) 


843 
36 


1400 
20 


4800 
6 


Si (mJy) 


125±10 


122±5 


69±5 


Reference 


Mathewson 
et al. 1984 


This 
work 


This 
work 



metal a bundances set to 0.5 solar as typical of the LMC; 
Russell & Dopitalll992l) . The statistical quality of the spectra 
was sufficient to fit one-component thermal plasma emission 
models. For the supernova remnant we used a single-temperature 
non-equilibrium i onisation colUsio nal plasma (NEI) model (in 
XSPEC, see Borkowski et al.ll200ll and references therein) with 
metal abundances fixed to 0.5 solar, yielding acceptable val- 
ues. Results from similar mod elling of X-ray spectra from S NRs 
in the SMC were published by van der Hevden etakl (|2004 and 
fFilipovic et al. (2008), allowing direct comparison. 

2.3. Tine MCELS opticai surveys of tine LMC 

The Magellanic Cloud Emission Line Survey (MCELS) was car- 
ried out from the 0.6 m University of Michigan/CTIO Curtis 
Schmidt telescope, equipped with a SITE 2048 x 2048 CCD, 
which gave a field of 1.35° at a scale of 2.4" pixel"'. Both 
the LMC and SMC were mapped in naiTow bands coiTespond- 
ing to Ha, [Oiii] (^=5007 A), and [Sii] (/l=6716,6731 A), plus 
matched red and green continuum bands that are used primarily 
to subtract most of the stars from the images to reveal the full 
extent of the faint diffuse emission. All the data have been flux- 
calibrated and assembled into mosaic images, a small section of 
which is sho wn in Fig. |5l Furth er details regarding the MCELS 
are given byi Smith et al.l (|200 6|). Winkler et al. (in preparation) 
and at http://www.ctio.noao.edu/mcels/ 

Additional observations were attempted on 25* September 
2009, using the 1.9-meter telescope and Cassegrain spectro- 
graph at the South African Astronomical Observatory (SAAO) 
in Sutherland. While the observing conditions were good, we did 
not detect any emission from SNR J0528-6714. 

3. Results 

3.1. Radio-continuum 

The remnant has a circular appearance centered at 
RA(J2000)=5''28'"18.5', DEC(J2000)=-67°13'49.2" with 
a diameter of 216"+5" (52.4+1.0 pc). This is reasonably 
consistent with its X-ray extent as indicated in Fig. |4] where the 
3cr radio-continuum contour is drawn on an X-ray image. The 
point-like X-ray source is located very close to the centre of the 
radio emission while the diffuse X-ray emission is brighter on 
the south-eastern side. 

A flux density measurement was made at 6 cm, resulting in a 
value of 69 ±5 mJy (lable[Ili. Using the 20 cm mosaic image de- 
scribed in Hughes et al.l (l2007l) . we made a new estimate of the 
flux density of SNR J0528-6714, with a value of 122 + 5 mJy. 
A spectral index (defined as 5 oc v") using the flux densities 
in Table [T] is estimated to g=-0. 36+0.09 (Fi g. |2]), which is in- 
dicative of older SNRs (iFilipovic et al.lll998h . We note that the 
36 cm and 6 cm observations are only interferometer measure- 
ments, whereas the 20 cm observations consist of both single 
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Fig. 1. ATCA observations of SNR J0528-6714 at 6 cm (4790 MHz) overlaid with fractional polarised intensity. The ellipse in the 
lower left corner represents the synthesised beam of 39"x31", and the line below the ellipse is a polarisation vector of 100%. The 
sidebar quantifies the pixel map and its units are Jy/beam. 



dish and interferometer measurements. This could lead to a flux 
density underestimate at these two wavelengths. Previous spec- 
tral index es timates were based s olely on single dish data, (a=- 
0.79±0. 15: iFilipovic et al.lfT998h . We note that this is a highly 
confused region, and that as our new e stimate uses higher res- 
olution data than iFihpovic et al.1 (Il998h . we have been able to 
exclude the confusing sources from the spectral index estimate. 

Linear polarization images at 6 cm were created using stokes 
Q and U parameters. This image reveals regions of moderate and 
somewhat irregular polarisation within the SNR. Without reli- 
able polarisation measurements at a second frequency we could 
not determine the Faraday rotation, and thus cannot deduce the 
magnetic field strength. 

The mean fractional polarisation at 6 cm was calculated us- 
ing flux density and polarisation: 

JSq 

P = 100% (1) 

Si 

where S q,S u and S / are integrated intensities for the Q, U and 
/ Stokes parameters. Our estimated peak value is P ~ 20% just 
off the centre of the remnant. Along the shell there is a pocket 
of uniform polarisation, at approximately 15% (Fig. [T]i, possi- 
bly indicating varied dynamics along the shell. This unordered 
polarisation is consistent with the appearance of an older SNR. 



3.2. X-Ray 

The EPIC spectra of the SNR (excluding the central source) are 
plotted in Fig. [3] together with the best-fit NEI model. The de- 
rived model parameters (LMC absorption, temperature kT and 
ionisation time scale t) are summarised in Table [2] together 
with inferred fluxes and luminosities. Flux and luminosity are 
given for the 0.2-2 keV band, determined from the EPIC-PN 
spectrum (which has the best statistics with a count rate of 
0.056 counts s"'). The intrinsic source luminosity with total Nh 
set to assumes a distance of 50 kpc to the LMC. The rela- 
tively large errors in the derived LMC absorption lead to large 
uncertainties in the (absorption corrected) luminosities for the 
soft X-ray spectrum of the SNR. Therefore, we give a luminos- 
ity range derived from fits with LMC Nh fixed at the lower and 
upper 90% confidence values. Similarly we derive a confidence 
range for the emission measure EM. 

Spectra extracted from the point-like source near the centre 
of the SNR yield low counting statistics and only a spectrum 
from EPIC-PN data with 5 bins between 0.3 and 2 keV could be 
obtained. The shape of the spectrum is similar to that of the SNR 
and we therefore used the same model for a fit to the PN spec- 
trum. The resulting best fit values for absorption column density 
and temperature are somewhat higher than derived for the SNR, 
but consistent within the statistical errors. It remains unclear if 
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Fig. 2. Radio-continuum spectrum of SNR J0528-6714. 
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Fig. 3. EPIC spectra of SNR J0528-6714 ([HP99] 498). The best 
fits using a single-temperature NEI model are plotted as his- 
tograms (black: FN, red: MOSl, green: MOS2). 





, . .. — 



5°29"20' ag^o' 40' 20' 2a"o" 
RA (J2000) 



Fig. 5. MCELS HQ'(top), [Sii] (Middle) and [Oiii] (bottom) im- 
ages of SNR J0528-6714 overlaid with an ATCA 6 cm 3o- 
(1.2 mJy) and lOcr (4.0 mJy) radio-continuum contours. 



the point-like source is caused by a (somewhat hotter) knot in 
the emission of the SNR or unrelated to it. No foreground star, 
which could have a similar X-ray spectrum, is found near the 
X-ray position. 

The spectrum of SNR J0528-6714 is characterised by a 
high absorption and a low temperature. Considering its large 
extent and low surface brightness, it is most likely an older 
remnant. The low tempera t ure is similar to 3 SMC SNRs 
(Ivan der Hevden et al.l 120041: iFilipovic et a l.' '2008), and cooler 
than a ny known brig h t LMC SNR (jW illiams et al. 2004; 
WilliamT&Ch3 12001 iBorkowski et"aLn2006l: iBambaetd.1 
2006t) . 



3.3. Optical 

There is no obvious association between any emission in the Ha, 
[Sii] or [Oiii] bands and SNR J0528-6714 (FigUll, confirming 
the earlier observation of Mathewson et al. (1985). The high- 
density group of stars KMHK 943 (Kontizas et al. 1990) is vis- 
ible on the eastern rim of the SNR J0528-6714. .Martavan et all 
(l2006h report several B and Be stars in this region. These 
other objects, and their interaction with the interstellar mat- 
ter (ISM), could explain the lack of an optical counterpart to 
SNR J0528-6714. Also, we cannot exclude the possibility of a 
dark cloud in the line of sight. 
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Fig. 4. XMM-Newton false colour image (red corresponds to X-ray intensity in the 0.2-1.0 keV band, green to 1.0-2.0 keV and 
blue to 2.0-4.5 keV) of SNR J0528-67 14 overlaid with an ATCA 6cm 3cr (1.2 mJy) and lOcr (4.0 mJy) radio-continuum contours. 



4. Discussion and Conclusions 

We have analysed our new and archival XMM-Newton, MCELS 
and ATCA observations in the direction to the SNR J0528-67 14. 
We conducted the highest resolution X-ray and radio-continuum 
observations to date of SNR J0528-6714. From radio observa- 
tions we found a diameter of 52.4±1. pc, at the top end of the 
size distribution dBadenes et al.ll20Tol) . a radio-continuum spec- 
tral index a=-0.36+ 0.09, which is indicative of an older SNR 
dFihpovic et al.ll998h . The unordered polarisation features, most 
likely due to a weakening of the remnants magnetic field, also in- 
dicate an older age for SNR J0528-6714. The almost perfectly 
circular appearance, along with the lack of radial syrametr y clas- 
sifies this object as a Type la SN event (iLopez et al.ll2009 *). 

However, the most unusual feature is the la ck of any opti- 
cal counterpart and very faint X-ray emission (iFilipovic et alj 
12008': 'Biliko va et alJl2007l and the references therein). While it 
is not unusual to not detec t an SN R in the optical, ( Filipovic et al. 
120011: IStupar et al.l |2004 l2005l) . the combination of faint X- 
ray and no optical detection is rare. This is in contrast to the 



M 33 findings of iLo ng et al.1 ([2010^, where no SNR is found 
with X-ray and radio emission, but not seen in optical. This 
may imply possible observational bias towards the optical tech- 
niques of SNR detection in external galaxies where the resolu- 
tion (and therefore s ensitivity) may pla y dominant role in SNR 
identification. Also, Long e t all (1201 Ol) found significant num- 
ber of M33 SNRs with similar (or lower) X-ray energies then 
SNR J0528-67 14 but their radio-continuum brightness is not ex- 
ceedingly high as in case of SNR J0528-6714. 

We detect regions of moderate and somewhat irregular po- 
larisation with maximum fractional polarisation of 15%. Also, 
we argue that as with the majority of other SNRs in the MCs, 
this intrig uing SNR is most likely in the adiabatic phase of its 
evolution ("Pa vne et alll2008l) . 

The nature of the point-like X-ray source near the center 
of the remnant remains unclear. It does not coincide with the 
maximum radio intensity which is found somewhat off-center 
to the East. It is also unclear if the radio maximum is truly 
a point source or just enhanced emission from the shell. This, 
and the lack of hard X-ray emission make the explanation that 
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SNRJ0528-6714 is a Pulsar Wind Nebula (PWN) very un- 
likely. Also, there is no reported radio pulsar in this region 
(Manchester et al. 20051). 
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